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1 Technical  /Hummary 

ABSTRACT  / 


An  algorithm  for  best  approximating  in  the  s\ip-norm  a function 


f e C[0,1]  by  functions  from  tensor-product  spaces  of  the  form 


IT.  9 C[0,1]  ® C[0,1]  8 IT  , is  considered.  For  the  case  k = i = 0 

K X> 

the  Diliberto  and  Straus  algorithm  is  knovm  to  converge.  A straight- 


forw£u:d  generalization  of  this  algorithm  to  general  k,£  is  formulated, 

and  an  example  is  constructed  demonstrating  that  this  algorithm 

2 2 

does  not  converge  for  k + £ > 0 . 

AMS(MOS)  Subject  Classifications:  41A50,  41A63 

Key  Words:  Algorithm,  Tensor-Product  Spaces, 

Best  Approximation  by  Polynomials 

Work  Unit  Number  6 - Spline  Functions  and  Approximation  Theory 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


Often  it  is  desireUt>le  to  approximate  a given  function  as  closely  as 

possible  by  a member  of  a class  of  functions  that  are  simpler  to  evaluate. 

For  a general  continuous  function  of  two  vauriables  f(x,y)  a best 

approximating  function  of  the  simpler  form  h(y}  + g(x}  can  be  computed 

by  the  algorithm  of  Diliberto  and  Straus.  Since  such  an  approximation 

can  be  quite  feu:  from  the  approximated  function,  a better  approximation 

of  the  form  J h.(y)x  + ^ 9.  My  is  considered.  One  way  to  try 

i=0  ^ j=0  ^ 

to  construct  such  an  approximation  is  to  generalize  the  Diliberto  and 
Straus  algorithm  to  this  more  general  setting.  The  generalized  algorithm 
is  simple  in  the  sense  that  only  one-dimensional  h>est  approximations  hiy 
polynomials  hiave  to  be  computed.  In  this  note  it  is  shown  by  a simple 
example,  that  this  "natural"  generalization  camnot  h)e  expected  to  con- 
verge, and  therefore  other  methods  should  be  developed  . 
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A STRAIGHTFORWARD  GENERALIZATION 


OF  DILIBERTO  and  STRAUS'  ALGORITHM  DOES  KJT  WORK 


Mira  Richter-Dyn 


The  algorltln  of  riliberto  and  Straus  for  approxiatating  a bivariate  function  by  a 


sum  of  univariate  ones  proposed  in  19S1  fl] , has  been  recently  investigated  in  several 


works  [2] , [3] , [4] , where  convergence  and  various  properties  of  the  algoritin  are 


studied 


is  of  the  following  form 


2n+l 


It  is  proved  in  [11,  [3],  [4]  that  limllf  11  = lnfllf-4ll>  although  the  rate  of  convergence 

n-x»  4cm 

might  be  extremely  slow  [2] . Algorithm  (2)  can  be  interpreted  as  a sequence  of  repeated 


f(x,y),  regarded  alternately  as  a function  of  x and  as  a function  of  y . More  specif 


ically,  let  J be  the  operator  associating  with  f(x,y)  c CtO,ll  the  function 


(J  f)  fy)  f.  C[0,1],  with  (J  i)  (y.)  the  constant  of  best  approximation  to  f(x,y  ) in 


the  sup-norm  on  10,1] > and  let  J be  defined  similarly  with  the  roles  of  x,y  inter' 


changed.  Then  (2)  can  be  rewritten  as 


*o“^'  ^2n+l"*2n"^x*2n*  ^2n+2*^2n+l'^y*2n+l 
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This  formulation  suggests  a straightforward  generalization  of  algorithm  (3),  namely 
best  approximating  f(x,y)  alternately  in  the  x and  y directions  by  polyncmials 
of  degree  k and  t respectively,  in  order  to  obtain  a best  approximation  to  f(x,y) 


(s.  denotes  the  space  of  all  univariate  polynomials  of  degree  ^ k.)  With  this  nota- 


algorithm  (3)  to  this  more  general  setting  is 


2n+l  2n  X 2n'  2n+2  2n+l 


2n+l 


where 


f) (x  ,y)  is  similarly  defined 


In  the  following  we  present  a simple  example  demonstrating  that  algorithm  (S)  for 


construct  a f'mction  f(x,y)  such  that  llfll  > inf  Ilf-^N,  while  the  functions  {f  } 

♦^”o.i 

generated  from  it  by  (S)  with  k>0,i’>l  satisfy  Ilf  II  • llfll  for  all  n • 


Consider  f(x,y)  CtO,l]  subject  to  the  following  conditions 


elsewhere  in  [0,1] 


As  can  be  easily  observed 


f satisfy  (6).  Thus  algoritln  (5)  with  k>0(twl 


genera'es  !•  sequence  {f  } of  functions  satisfying  (6)  whenever  f.  satisfies  (6) 


'2' 


In  order  to  verify  that  Ilf II  > inf  llf-^ll,  it  ie  sufficient  to  show  that  there 

♦‘"o.i  2 

does  not  exist  a bounded  linear  functional  li  e (ClOtl]  ) ' « w 0,  such  that 


Indeed  any  v¥0  with  property  (8)  is  necessarily  of  the  form 


with  r > 0 


namely  a linear  combination  of  function  values  at  extremal  points  of  f . Moreover 
condition  (7)  implies  that  u can  be  written  as  a linear  combination  of  first  differences 
in  the  X direction  so  as  to  vanish  on  all  functions  of  the  form  h(y) , and  as  a linear 
conbination  of  second  order  divided  differences  in  the  y direction , so  as  to  vanish 


on  all  functions  of  the  form 


These  characteristics  of  u are  consistent  with  the  special  structure  of  the  15 


extrmnal  points  of  f # as  given  in  (6)  , only  if  r=“14  in  (9) . Then  v can  be  written 


denotes  the  second  order  divided  difference  of 


linear  combination  of 


points  of  f with  x* 


first  differences  in  the  x direction  only  if  c^ 


which  admits  only  the  trivial  solution 
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